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Abstract 

The use of dyadic interaction between agents, in combination with homophily (the principle that 
"likes attract") in the Axelrod model for the study of cultural dissemination has two important 
problems: the prediction of monoculture in large societies and an extremely narrow window of 
noise levels in which diversity with local convergence is obtained. Recently, the inclusion of social 
influence has proven to overcome them (A. Flache and M. W. Macey, arXiv:0808.2710). Here we 
extend the Axelrod model with social influence interaction for the study of Mass Media effects 
through the inclusion of a super-agent which acts over the whole system and has nonnull overlap 
with each agent of the society. The dependence with different parameters as the initial social 
diversity, size effects, Mass Media strength and noise is outlined. Our results might be relevant 
in several socio-economic contexts and for the study of the emergence of collective behavior in 
complex social systems. 



I. INTRODUCTION. 



Social human systems are constituted by human agents which interact with a limited 
number of peers. This perhaps apparently random procedure yields to global regularities: 
there are order- disorder transitions (spontaneous formation of cultures, the emergence of 
consensus about different issues, spontaneous formation of a common language), examples of 
scaling and universality, to mention just a few. These phenomena are now been under intense 
studies through different techniques. One of them is the so called agent-based modelling 

il-bi- 

The interest in the application of such techniques has grown rapidly mainly as a result 
of the increasing availability of computational capabilities. Today, computer simulations 
have become an excellent way to model and understand social processes. The usefulness of 
social simulation modeling results as much from the process (problem specification, model 
development, and model evaluation) as the product (the final model and simulations of social 
system dynamics). Nevertheless, social simulations need to be a theory-guided enterprise. 
Its results will often be the development of explanations, rather than the prediction of specific 
outcomes [lj]. 

Agent-based simulations allows us the possibility to develop a new way of thinking, par- 
allel to analytical procedures, based on the idea about the emergence of complex behaviors 
from relatively simple activities [lj]. It has great potential to assist in discovery of simple 
social effects when social scientists build very simple models that focus on some small aspect 
of the social world and discover the consequences of their theories in the "artificial society" 
they have built. 

One of those simple models already built is the well known Axelrod model for the study of 
dissemination of culture jl-sj]. He argued that culture "is something people learn from each 
other", and hence something that evolves through social influence. To study the process of 
cultural diffusion Axelrod built a model based on two simple assumptions: 

1. people are more likely to interact with others who share many of their cultural at- 
tributes, and 

2. these interactions tend to increase the number of cultural attributes they share (thus 
making them more likely to interact again). 
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The Axelrod model has been exhaustively implemented to study a great variety of prob- 
ems: the nonequilibrium phase transition between monocultural and multicultural states 



10J, the cultural drift driven by nois e [ill . Il2j. nominal and metric features jl^, propa- 



ganda 



15], time evolution dynamics |l6lj. the resistance of a society to the spread of foreign 



cultural traits [17|, finite size effects [18[, the impact of the evolution of the network structure 



with cultural interaction [ljj, and mobility of social agents [20|, among others. 

The dynamical influence on the society is described in the original Axelrod model using 
a dyad interaction which, at each time step, only comprises the source and the target of 
influence in a particular interaction. As discussed by A. Flache and M. W. Macy in Ref. 
2l| . the influence a person feels from the society is a social phenomenon that cannot be 
reduced to the interactions within a dyad given by source-target couple of persons because 
the social pressure on the target to adopt an opinion is proportional to the number of people 
that the target perceives are supporting this opinion. On the contrary, the Axelrod model 
assumes that influence is interpersonal (dyadic). We are then in front of a two completely 
different assumption with respect to the way people interact in society: 

1. dyadic interaction where it is supposed that two people in a relationship interact in 
isolation from others and 

2. social influence which is multilateral and involves all network neighbors simultaneously. 



In Ref. [2l| the authors explained that the central implications of Axelrod's model 
profoundly change if the dyadic interaction is changed to social influence. The authors have 
shown that the combination of social influence with homophily (the principle that "likes 
attract") solves two important problems which contradict the empirical pattern: 

1. the original Axelrod model predicts cultural diversity in very small societies, but mono- 



culture in larger ones 



30j and 



2. when cultural perturbation is present, diversity is obtained only in a very narrow 
window of noise level and this window decreases with increasing population size 
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On the other hand, one of the extensions of the Axelrod model that have been done 



consists of the study of Mass Media effects over the society [22|-|25| . Surprisingly, the Mass 
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Media has been found to be a destabilizing 



state 



261 ] . The exception are the works 



actor which drives the system to a multicultural 



27-291 where different models has been introduced. 



The aim of this work is to investigate the effects of Mass Media over the society when 



the social influence (in the sense of Ref. 



21] ) is considered. This paper is organized as 



follows: in section [TT] a brief description of the Axelrod model and its results are explained, 
in section II I II it is introduced the model which considers social influence, in combination 
with homophily and Mass Media Our results are analyzed in section IIVI and finally some 
conclusions are outlined. 



II. A BRIEF DESCRIPTION OF THE AXELROD MODEL. 

The Axelrod model consists of a population of agents, each one occupying a single node of 
a square network of size L and area Lx L. The culture of an agent is described by a vector of 
F integer variables {erf} called features (f = 1, F). Each feature can assume q different 
values between and q — 1. These are the possible traits allowed per feature. In the original 
Axelrod model the interaction topology is regular bounded (non-toroidal). Each agent can 
interact only with its four neighbors (a von Neumann neighborhood) which are the most 
closer (only one step distant from the target agent of influence) without crossing the borders. 
Initially, individuals are assigned a random culture and the parameter q, which defines the 
possible traits in each cultural dimension, can be seen as a measure of the initial disorder 
or cultural variety in the system. In the temporal dynamics of the model at each time step 
t, the cultural profile of a random selected agent i may be updated through the interaction 
with a randomly chosen neighbor j. The probability of this interaction is proportional to 
the corresponding overlap of their cultural profiles (the amount of features with identical 
traits). When interacting, agent j influences agent % causing the last to adopt j's trait on 
a feature randomly chosen from those that they do not share. The process outlined above 
continues until no cultural change can occur. This happens when every pair of neighboring 
agents have cultures that are either identical or completely different. 

In the model it is explored the aggregate behavior by studying the spatial distribution 
of the emergent cultural regions: sets of spatially contiguous agents who share an identical 
vector of culture. Two final possible states are possible: only one cultural region is obtained 
or multiple cultures are obtained separated by a boundary. These states are called mono- 
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cultural and multicultural (or global, polarized), respectively. A first order phase transition 
is obtained for F > 2 between the monocultural to a multicultural state for increasing value 
of the initial diversity q 9|. 

The Axelrod model has shown how diversity with local convergence can be obtained. 
Nevertheless, as mentioned before, there are two results of the model which contradict 
empirical patterns: (1) the Axelrod model predicts cultural diversity in very small societies, 
but monoculture in larger ones and (2) when cultural perturbation is present, diversity 
is obtained only in a very narrow window of noise level and this window decreases with 
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ll|. This has been justified as a consequence of the dyadic 



increasing population size 
interaction considered between the agents of the society [21]. 



III. THE MODEL FOR SELF-INCLUDED SOCIAL INFLUENCE WITH MASS 
MEDIA. 



Following the proposition made by A. Flache et al in Ref. [2l|] we have developed a 
new interaction procedure in the Axelrod model to study the effects of the social influence 
between agents. Our approach is similar, but not exactly equal, to that proposed in Ref. 



211 ] . In our case, at a given time t an agent i is randomly chosen for possible influence. 
Around any agent, we consider a neighborhood defined by a parameter p, which controls 
how many agents can be visited (in unitary steps) in each direction starting from that agent 
until the perimeter of the neighborhood is reached. One unit-step allows to move only in 
north-south or east- west direction. When p—1 the neighborhood obtained is therefore the 
well-known von Newmann neighborhood N±, which only includes the four closer neighbors 
of an agent (Aq = {31,32,33,34})- In our case, p can take higher values to define bigger 
neighborhoods. In general, for a given value of p it is obtained N p = {ji,---,j a } where 
a = a(p) = 2p(p + 1). This give four neighbors (a = 4) for one step neighborhood (p — 1), 
twelve neighbors (a = 12) for two step neighborhood (p = 2) and so forth. 

The social influence is defined as follows: when an agent i is selected for possible influence, 
all of its neighbors (according to the value of p) are included in the set of influence Ij with 
probability pij a = Oij a /F, where Oij a is the cultural overlap between agent i and j a . The 
set Ij obtained by this procedure becomes a set of traits' influence over agent i's traits. In 
our model we have decided to include also agent % into the set I, with probability pa = 1. 
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FIG. 1: Social interaction between agent i and its neighbors. Each neighbor is included with 
probability pij a into the set of influence Ij. The agent i is also included with probability pa = 1. 
The agent s, included with probability pi s according to Eq. (UJ represents a super-agent which 
acts over the whole society. 

This allows a self-reflection of agent i about its own traits when comparing with the traits 
of its neighbors. Then, each agent is aware of its own traits when being influenced by its 
neighbors and its own traits count when deciding to change, or not, its corresponding value 
for a different one. 

Figure [1] illustrates the social influence procedure. In the example shown, the agent i 
selected has cultural traits (1,1,1) and the neighborhood is formed by the one-step neighbors 
(von Neumann neighborhood). We have also supposed that all the four neighbors have 
entered into the set of influence Ij. In general, neighbors are included according to the 
probabilities Pij, computed for each neighbor. The last agent to be included on the set Ij is 
an agent shared by all the agents on the network. Without loss of generality this super-agent 
s is set up with culture zero in all of its trait values and never change them. This global 
agent represents a Mass Media influence which acts over the whole society. Its inclusion in 



this model follows the same rules explained on Ref. 



29[ , according to which the probability 



of interaction between this super-agent with an agent i of the society is given by 

O is + e 

p« = -WT7 ( } 

where Oj S is the overlap between the agent i and the super-agent s. The parameter e(> 0) 
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can take real values and it is a measure of the Mass Media strength to act over all agents in 
the society. Then, the super-agent s has always a non-zero probability (min(pj s ) = e/(F + e)) 
to be included in the set of influence of the agent i. This probability is proportional to the 
strength parameter e and its non-zero value means a clever Mass Media which designs its 
publicity actions to have always something in common with the target of influence. 

Once the agent % is randomly selected and the current influence set Ij is established 
(including i and s), an agent z's feature is selected at random with a trait value different 
from the corresponding value of at least one agent of the set Ij in the same feature. The 
decision of which trait value agent % adopts (if any change takes place) is done in such a way 
that agent i imitates (or copy) the most common trait in the set Ij for the feature selected. 
If, for example, in the situation depicted in Fig. [JJ the feature selected on agent i were the 
second (that in the middle), this trait value (0^2 = 1) will be changed to 0^2 = due to 
highest frequency of this last value in the set Ij. If the feature selected were the third (from 
top to down), the value 0^3 = 1 will remain due to the highest frequency of this value on the 
set Ij. The last possible case represented in Fig. [1] occurs if the first feature were selected. 
In this case, on the set of influence Ij, we have the same frequency of appearance for the 
values and 1. In this situation, one of the two values is chosen at random and then, the 
trait value 0^1 = 1 in agent % could remain with probability 1/2 or changed to 0^1 = with 
the same probability. Then, in general the trait value of agent i is changed to the value of 
the corresponding trait with the highest frequency of appearance on Ij. If there are more 
than one, it is changed to one of them with equal probability. The trait value of agent % 
does not change if this is the only one with the highest frequency of appearance on Ij. 

At the same time, we have also included interaction errors as well as copying errors, just 



as in Ref. 



2l|. The interaction error relaxes the previous deterministic procedure used in 



former works of Axelrod model when deciding the possible interaction between two agents 
(dyadic interaction). Both copying and interaction errors can randomly alter the outcome 
of an interaction event. The interaction error acts over the selection procedure as follow: 
if the normal procedure of selecting an agent j for being included in the set of influence Ij 
results in its inclusion, then with probability r' the neighbor j is removed from Ij. If the 
neighbor has not been selected into the influence set, then with probability r' the agent 
j will be included into Ij. This error creates the possibility that cultural influence occurs 
across the boundary of two disconnected cultural regions if a neighbor with zero overlap is 
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included in the influence set. This error can also reduce the social pressure against adopting 
different trait values when a cultural identical neighbor is excluded from the set of influence 
and increases the possibility for the target agent % to adopt trait values from a completely 
different culture. On the other hand, the copying error acts after agent % has adopted (or 
not) a different trait value and the new trait value has been already set up. In this case, 
the corresponding trait value adopted by agent % is changed to a new randomly selected 
value with probability r. Notice that the new value randomly generated could be that one 
the agent i already has. In general, both interaction and copying errors are conceptually 
different, but for simplicity we have used here the same value r = r'. Additionally we have 
implemented periodic boundary conditions (toroidal society) to avoid boundary effects. 

In the case of dyadic interaction and in the absence of noise, when analyzing the possible 
final absorbing states, it is simpler to establish it by checking when each agent has full or null 
overlap with each one of its neighbors. In the case of social influence the problem is more 
involved. However, it is also possible to establish some technical conditions to be checked 
to see if an agent is active (i.e., can interact with its neighbors according to the dynamical 
rules established in the model) or not taking in consideration the set of influence L> When 
all the agents of the society are inactive, then an absorbing state is obtained. 

In the present study we have also implemented the procedure developed in Ref. 30 1 
where a list of active agents is built. In this way, instead of randomly selecting agents of 
the society in each time step, the agents are randomly chosen directly from this list. This 
procedure strongly increase the efficiency of the dynamical evolution of the system and it 
allows to save computational time. Therefore, when the system is initialized by randomly 
assigning different cultures to each agent of the society, the first list of active agents is 
built. Next, at each time step, when the influence is established and an agent of the society 
changes its cultural value, the list of active agents is updated analyzing the agent itself and 
all its neighbors to check which of them are now active. The dynamical iteration keeps on 
until the length of the list of active agents is reduced to zero. It is worth noticing that in 
our case (where social influence is established), some runs did not settle into a final, well- 
defined absorbing state. In these cases the list of active agents reduces to one element and 
each time this agent changes its cultural values and becomes inactive, one of its neighbors 
becomes active. This propagation goes on indefinitely. We have neglected these cases from 
our calculations and we have only considered runs which finish in a precise well-defined final 
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absorbing state. Moreover, when noise is included, it is not possible in general to define a 
final absorbing state and the dynamical evolution of the system is stopped by defining some 
criterion related with a definition for the stationary state the system reaches. In our case, 
we have included noise only in those agents which are active on each time step, and do not 
in the rest of the society. This has allowed us to reach final absorbing states even when noise 
is present. 

IV. RESULTS 

We have implemented different computational experiments to study the Mass Media 
effects in the Axelrod model when the social influence is the mechanism at work for agents 
interaction. As there is no qualitative difference on the dynamical behavior of the system 
for F > 2, we have set F = 3 in our study, a network size with L = 50 (2500 agents 
unless otherwise stated) and averages have been taken over 200 different initial random 
configurations. 

In Fig. [2] we represent the normalized value of S (the amount of agents in the biggest 
culture domain when the system reaches the final absorbing state), as a function of the 
initial diversity q with no noise (r = 0). Panel a) corresponds to calculations with one-step 
neighborhood (p = 1) while in panel b) the calculation was done for p — 2. This means that 
the amount of neighbors of each agent in panel a) are only four plus the super-agent while in 
case b) each agent has now twelve neighbors beside the super- agent. Different values of the 
super-agent strength e are included. It is also included in dotted-line the parameter 8 which 
is the amount of information of the super-agent on the society. This parameter is obtained 
counting along all of agents and their corresponding traits if the information is equal to that 
of the super-agent in the same feature. The final count is divided by F to scale the value in 
a way that can be shown in the same scale as the parameter S. 

It can be seen that for low values of q the system reaches a close-mono cultural state for 
any value of the strength e. These are rather trivial cases because these extremely low values 
of q mean a very low initial diversity of the system and a final close-monoculture state is 
then expected. 

A monoculture state is also induced for high enough values of q but the system is now 
more sensitive to the values of e than in the region of low q values. In this case, higher 
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FIG. 2: Final absorbing state in the Axelrod model with social influence as a function of q for 
different values of e. In dots-symbol line is the fraction of the biggest culture while in dotted line 
is the fraction of Mass Media information on the society (only in panel a)). The parameter A is 
the area of the system (A = L 2 ). a) Calculation done including one-step neighbors (von Neumann 
neighborhood), b) Calculation done including two-step neighbors. The calculation is done with F 
= 3, L = 50 and averages are taken over 200 different random initial conditions. 

values of e induce a stronger multicultural final state given by higher values of < S >. To 
explain this result, we note that higher values of q mean there are initially a higher degree of 
cultural diversity on the society and this is reflected in sets of influence I with low number 
of neighbors, i.e., each set of influence I, will be frequently composed by the own agent i 
and by the super-agent s, and with low probability by the neighbors since they probably do 
not share any of their trait values with agent i. The probability of the super- agent to be 
included in the set of influence I increases for increasing value of e. For higher q values, in 
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cases where Ij is formed by agent i and super-agent s, the last one will be able to introduce 
its own value on agent % with probability 0.5. The iteration of this process in time drives 
the system to a close monocultural state. 

A more interesting situation occurs for middle values of q where it is observed a minimum 
of the < S > I A values as a function of q, more pronounced for e = 0.05 and q = 10 in Fig. [2] 
a) and e = 1.00 and q — 15 in Fig. |2jb). For values of q close to the corresponding minimum, 
the initial diversity is such that besides agent i and s, some of the agent's neighbors are also 
included in the set of influence I. Then the interacting dynamics involves higher cultural 
diversity on the set of influence Ij and the super-agent fails to induce a strong monocultural 
final state. 

It is interesting to note that when comparing Fig. [2] a) and b) it can be seen that the 
increase of the number of neighbors included in the social interaction decreases the size of 
the biggest culture in the absorbing state and then a more pronounced multicultural state is 
reached for the same values of q and e. The value of q where the minimum of S is attained 
also increases. This is a consequence of a direct competition between the higher diversity 
on the set Ij and the homogenized influence of the super-agent. In Fig. [2jb), the amount 
of neighbors to be analyzed for inclusion on the set Ij is twelve. Three times the case in 
Fig. |2j a), which is only four. Then with p — 2, at any value of q there will be, with 
higher probability, bigger diversity of trait values than in the case of p — 1 and, therefore, 
a weakening of the homogenizing effects of the super-agent is expected. This decreases the 
possibility of the super-agent to drive the system to a monocultural absorbing state because 
the probability for a trait value of the super-agent to appear with the highest frequency on 
the set of influence is lower and lower values of S are obtained. Then, when considering more 
neighbors in the social interaction the higher local diversity reinforces the final multicultural 
state, even with the presence of a mass media. 

On the other hand, when comparing the dots line for < S > with the dotter line for 
< 6 > in Fig. |2j a) it can be seen that they both coincide for all values of q when e has 
higher values (e = 1.00 and e = 3.00). It means that the information of the super-agent is 
present only in the biggest culture of the society. Nevertheless, for small values of e, the 
two curves coincide for higher values of q, but the values of < 5 > are higher than the 
values of < S > for middle values of q. Then, in this situation although the final state does 
not correspond to a strong monoculture, the information of the super-agent on the society 
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FIG. 3: Final absorbing state in the Axelrod model with social influence as function of the Mass 
Media strength e for different values of q. a) Calculation done including one-step neighbors (von 
Neumann neighborhood), b) Calculation done including two-step neighbors. The parameters used 
for these calculations are the same as in Fig. 

include the bigger culture and also agents which belong to smaller culture sizes. 

Figure [3] shows a second computational experiment. Here, we have included the values 
of the same parameter < S > /A, but now as a function of the super- agent strength e 
for three different values of q. In this case, noise has also been neglected. One and two- 
steps neighborhoods are included in panels a) and b) respectively. For the three values of 
q shown, a monoculture is strongly induced as e increases due to the constant presence of 
the super-agent on the set of influence. The order of the curves according to the value 
of q is in correspondence with the order of the curves in Fig. |2] and the values of e. As 
e increases the super-agent becomes a factor of "normalization" making constant its own 
values in time along the dynamical evolution of the system. Even for q = 10, when the set 
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of influence includes with high probability some of the neighbors of the target agent and the 
super-agent traits values are not on the majority frequency of appearance, the latter succeed 
in inducing a monoculture as the strength e increases. In particular, for e > 1, an almost 
monocultural state with < S > /A > 0.7 is already induced, i.e., at least 70% of the society 
belongs to the biggest culture for all values of q included. For e = 1, we have min(p is ) = 
0.25 (this is obtained considering null overlap between the super-agent s and an agent i) but 
the super-agent can induce the 70% of the society to belongs to the bigger culture with its 
own information (according to the value of the parameter 5 shown in Fig. [2]), for any initial 
diversity q. The maximum value of e included on Fig. [3] is e = 27, which gives min(pj S ) 
= 0.90. Furthermore, when comparing in Fig. [3] the panels a) and b) the same results are 
obtained as when comparing in Fig. [2] a) and b). In this case the inclusion of a two-step 
neighbors in the social influence dynamics of the system strongly decreases the value of S 
for the case q = 10 and 30. 

In order to get additional insights on the system dynamics, we next consider the same 
model but including noise. In Fig. H]it is reported the dependence of the parameter < S > /A 
with respect to different noise levels for several values of q, e and L. In panels a) and b) we 
have considered one- and two-steps neighborhoods, respectively. The noise included ranges 
between and 0.45. As can be seen, social influence makes absorbing states more stable to 
a bigger range of noise level than dyadic interaction, as already reported in Ref. 2l|. In 
general there is no qualitative change for at least three order of magnitude for all the values 
of q, e and L used. For higher strength of Mass Media (e = 3.00, mm(pi s ) = 0.5) the final 
absorbing state remains almost monocultural for four orders of noise level and finally for 
noise values higher than 0.1 the system is driven to a full monocultural state reached at r 
= 0.45 approximately. The situation is different for low values of super-agent strength (e = 
0.05). In this case the final absorbing state is also stable to different noise levels, but only 
for three orders of magnitude. In panel a) of the Figure it can be seen that for q = 5 and 
noise level higher than r = 0.01, the system is driven to a monoculture. For increasing noise 
level at q = 10 and 30 (for e = 0.05) there is first a reinforcement of the multicultural state, 
given by decreasing values of < S > /A (stronger induced for q = 30), and later for higher 
values of r it is induced a monoculture. All the calculations shown in black circles have been 
done with a network size of L = 50. In order to study if the stability of the final absorbing 
state to noise is robust with respect to the network size, we have also explored the case L 
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FIG. 4: Dependence of the final absorbing state with respect to different noise levels and different 
values of q. In solid line it is included the case with e = 3.00 while dashed lines are the cases with 
e = 0.05. In black circles are reported the final absorbing states calculated with a network size 
of L = 50, while in white triangles the case with L = 10. a) Calculation done including one-step 
neighbors (von Neumann neighborhood), b) Calculation done including two-step neighbors. 

= 10. Quantitative differences arise only for q = 10 and 30, for e = 0.05 (white triangles 
in Fig. H]a)). In general it is observed that the final absorbing state is stable to the same 
range of noise independent of the network size. The difference obtained, consisting in the 
higher values of < S > /A for q = 10 and the higher minimum for q = 30 at r = 0.05 (both 
for e = 0.05) are due to the lower size of the network (L = 10). There are no quantitative 
differences for the other parameters due to the rather trivial case of q = 5 (for any value of 
e > 0) and for the high value of e used in the solid lines. 

In Fig. H]b) we have included the same results of panel a) calculated with L = 50, but 
now considering two-step neighbors. The results are qualitatively the same for e = 3.00 
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(no matter the values of q) and for q = 5 and e = 0.05. For the case of e = 0.05 and q = 
10 and 30 the system ends in a multicultural state for almost all the range of noise levels. 
Nevertheless, for sufficiently high noise level the system is abruptly driven to a monocultural 
state, as in all the others cases. 

Here we advance an analysis for possible explanation of the results obtained on this 
experiment. There are three parameters involved: the initial cultural diversity q, the strength 
of the super-agent e and the noise level r, and two different error mechanisms: the copying 
and including error. For low values of e the super-agent will appear on the set of influence I» 
with low frequency. This frequency increases for increasing e. For low values of q the majority 
of the neighbors of the target agents will be included on I, because of the low cultural 
diversity, while the opposite occurs for large enough values of q. Nevertheless, as can be seen 
from the Figure, this complex interplay does not has an important impact on the dynamics of 
the system along a wide range of noise levels and the value of < S > remains almost constant. 
When r > 0.1 the copying error becomes the dominant mechanism because independently 
of the cultural trait an agent has copied from its neighbors, the copying error changes it 
to any one randomly selected and as it was already said, it deletes the boundaries between 
different cultural regions. It is also of importance that the error mechanism stops when the 
agent becomes inactive. Hence, the social influence is the mechanism which allows an agent 
to be active/inactive and it is also the mechanism which switch on/off the copying error. 
Then, as the copying error connects two completely different cultures the homogenization is 
favored and agents become inactive (and noise stops) when its culture is completely equal 
to its neighbors. Therefore the final monoculture obtained at r = 0.45. 

At lower values of noise a non monotonic behavior is obtained for q = 10 and 30 when 
the strength e has low values (e = 0.05 in this case). The value of < S > first decreases for 
increasing noise reinforcing the multicultural state. This effect is much pronounced at q — 
30. In these cases, the initial diversity makes that the set Ij be formed only by the target 
agent and some of its neighbors. More neighbors will be present on the set of influence for 
q — 10 than for q = 30. For increasing noise the including error makes the set Ij populated 
by both agents from cultures which share some traits with the target agent and also with 
agent neighbors which do not. The super- agent will also be included. This interplay seems 
to produce strong local convergence and it drives the system to a multicultural state for a 
range of noise between 0.01 and 0.1. Finally, as already noted, the copying error drives the 
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Area 



FIG. 5: Culture-area relationship for q = 10 and different values of e. In circles are included the 
calculations for the corresponding value of e indicated with labels. In square are represented the 
calculations done according to the relation e = k(L — 5), with the value of k indicated with a 
label. F = 3, L = 50 and 200 different random initial conditions were used. The area is obtained 
according to A = L 2 . 

system to a monoculture state for even higher values of r. 

Finally, we study the influence of the Mass Media in the culture-area relation when the 
interaction on the society is social. An extensive examination of the culture-area relation 



on the Axelrod model with dyadic interaction was done in Ref. [30|. The authors obtained 
a non-monotonic behavior for the culture-area relation for q values below the critical value 
where the transition occurs (q < q c ) while for q > q c the number of cultures G first increases 
in a power-law dependence as G ~ A x with x = 1 and then gradually flattens when the area 
becomes of the order of the maximum number of cultures q F . In the limit L — > oo there 
are only two possible outcomes: for q < q c a single culture dominates in an ordered regime, 
while G — > q F and all the cultures are represented in the network in a full disordered regime. 
The authors proved the transition between these two regimes to be discontinuous because 



when the model is changed including instead 



G jumps from 1 to q at q = q c . Contrari 
of dyadic interaction social influence, Ref. 
which produces a multicultural final state only for small societies 30] is overcome and only 



21] , the problematic result of the Axelrod model 
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multicultural final states are obtained. In case the Mass Media is included according to the 



model used in Refs. 



23 



24J and 25] which take the interaction between agents as dyadic, it 



was proven in Ref. [26[ that even a vanishingly small value of the Mass Media probability to 
interact with agents is sufficient to destabilize the homogeneous regime for very large lattice 
sizes, contrary to the usual commonsense. 

Here we have calculated the amount of different cultures present in the final absorbing 



state indicated by the parameter G. We have followed the definition of Ref. [30( where 
all the different culture configurations obtained are counted without paying attention to 
simple connected regions. As there are q F different possible cultures, the parameter G/q F 
is normalized and its maximum value {G/q F = 1) means a completely multicultural state 
while low values are related with monocultural final states. 

In Fig. [5] it is included with circles the final cultural diversity G for a fixed value of the 
Mass Media strength e. It can be seen that at e = the amount of cultures obtained in 
the final absorbing state increases for increasing value of the area and reaches the maximum 
multicultural state for sufficiently high value of A. The calculation done with any other value 
of q gives the same dependence of the area which includes the same slope and only a parallel 
shift. The other curves were then not included for simplicity. We have then qualitatively 
reproduced the results obtained in Ref. [21|. In our case, a calculation of the slope yields 



x = 0.42 ±0.02, which is different to that reported in Ref. [30j. We think that the deviation 
is due to the difference that arises when dyadic interaction or social influence is included. 

When the super-agent is included the values of G are lower for higher strength e at the 
same value of the area, but the slope of the curve is the same as when the super-agent is 
not present. The super- agent also prevents the system to reach the full multicultural state 
(< G > jq F = 1) for higher values of the area. This can be seen on the curve for e = 0.05, 
which saturates at values of A > 10 4 . The same seems to occur with e = 3.00. Calculations 
taking in consideration two-step neighbors were also carried out with the same qualitative 
results (including slope and saturation). Only a shift to decreasing values of G were obtained 
for the case of e = 0.00 and 0.05, but the shift finally disappears for e = 3.00 meaning that 
for strong enough Mass Media the relative importance of the amount of neighbors including 
in the set of influence is weak (at least for the value of q = 10). 

Furthermore, on the dynamics of the system there are two parameters which compete to 
produce opposite effects for increasing values of each one. For increasing area of the network 
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the system tends to reach a multicultural state while an increasing value of the super-agent 
strength e pushes the system to a more cultural homogeneous absorbing state. In order to 
study the relative weight these parameters have over the system, we have also calculated 
the final absorbing state when both area A and strength e increases. To accomplish this 
purpose we have established a dependence of the strength e to the area given by the following 
relation: 



In Fig. [5] it is included, with square symbols, the results for k = 1 and 10. The values 
of the diversity G obtained are higher for k — 1, than for k = 10, due to the lower rate 
of increase of the strength e and consequently a more multicultural final state is induced. 
Nevertheless, both curves increase for increasing area, which means that the area of the 
network has more weight than the strength e on the dynamics of the system, as expected if 
we examine Eq. (JTJ) . In our case, the maximum overlap between the super-agent and any 
agent of the network is max(Oj S ) = F — 3. Then, for e/F ^> 1 (which occurs rapidly for 
increasing value of A in Eq. EJ) it is obtained than pi S — > 1 and the relative increase of the 
probability Pi S is cancelled out with the effects of the increasing area. This explains also 
why the curves for k = 1 and 10 tend to the same values as the area is increased. 

Then, when social interaction is present according to the present model, increasing net- 
work size always drives the system to a multicultural state, while increasing super-agent 
strength prevents the system to reach the maximum possible of cultural configurations. Ad- 
ditionally, the saturation value of G seems independent of the value of the strength of the 
mass media, as far as our calculations have shown. 

We have also included in dotter line in Fig. [5] the analytical expression reported in Ref. 



which is the average number of cultures in the totally disordered configuration were A = L 2 
agents are randomly assigned with one of the q F available cultures. The expression is valid 
for A and q F very large. As can be seen, the prediction from Eq. (j3J overestimates the 
values for the cultural diversity as a function of the area. The results are more different the 
higher the values of the super-agent strength e are, since it is a factor that decreases the 
cultural differences between neighboring agents. 



e = k(A- 5) 



(2) 




(3) 
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V. CONCLUSIONS. 



Summing up, we have numerically studied an Axelrod model in which two different 
mechanisms are at work. The first takes into account that agents are prone to change 
their state according to a sort of social influence given by a (variable) neighborhood made 
up of their nearest (both in space and in cultural similarity) neighbors. The dynamics of 
the agents is also ruled by a Mass Media effect, represented in this case by an agent that 
keeps fixed its cultural traits all the time. The dynamics of the system is such that at 
low and high values of initial social diversity q, a monocultural state is attained, with a 
stronger dependency on the mass media strength (e) at large q values. The reason behind 
the formation of these monocultural states is the varying size of the set of neighbors with 
social influence: the higher the diversity of initial cultural traits, the less is the number 
of agents in the neighborhood of an agent. Admittedly, the more interesting behavior is 
for intermediate values of q, where the system dynamics attains a minimum in the size of 
the biggest monocultural cluster. At those values of q the initial diversity is such that the 
mass media has to compete with a larger number of neighbors of the agent and thus its 
information is not necessary on the majority. Then, in this case the mass media fails to 
drive the system to an homogeneous cultural state. Nevertheless, the increase of the mass 
media strength for a fixed network size always reinforces the monocultural state. 

On the contrary, the increase of the network size was proven to be a parameter which 
strongly drives the system to a polarized society where all the possible cultural configurations 
are present. The amount of final number of cultures also follows a power-law dependence 
when the social influence is at work, but the exponent value found here is lower than that 
reported for the case of the original Axelrod model with dyadic interaction. When mass 
media is present it only succeeds in preventing the full multicultural state for increasing 
size of the society and the amount of cultures obtained saturates. Moreover, the saturation 
value was found to be independent of the mass media strength in our calculations. 

We have also studied the model with noise representing errors in copying traits and in 
the formation of the social influence interactions. In radical difference with previous works 
where the noise is present on the whole society, we have included here noise effects only in 
those agents which are active according to the rules of the social interaction. This allowed 
us to reach well-defined final absorbing states and correspondingly a higher precision in the 
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description of the noise effects. The results show that social influence makes the system 
dynamics more stable against the presence of noise and that the latter only has a marginal 
influence on the general qualitative picture obtained without any errors. It is however 
worth stressing that noise has in general a positive effect in the formation of monocultural 
final states, giving rise to polarized societies for large enough values of it, independently 
of the other parameters ruling the size of the social influence and the strength of the mass 
media effect. Nevertheless, the fact that we have not distinguished between the two error 
mechanisms in this study does not allow to better understand the specific role each one 
plays on the dynamics of the system and further studies are needed where these two types 
of errors were considered as independents. 

We also think about the need to conduct more detailed studies exploring the implications 
of the social influence with respect to the dyadic one. The introduction of different models 
which better study the inter-relations between the dyadic interaction and social influence 
on a social system could be of importance to further elucidate the robustness of the results 
found here. This work is currently in progress. 

Acknowledgements 

A. H. R. thanks to the International Center for Theoretical Physics (ICTP), were the 
major part of this work was carried out. Y. M. is supported by Spanish MICINN through 
projects FIS2008-01240 and FIS2009-13364-C02-01. 



[1] N. Gilbert and K. G. Troitzsch, Simulation for the social scientist (Open University Press, 
England, 2005). 

[2] Artificial Societies, edited by N. Gilbert and R. Conte (Routledge, London, 1995). 
[3] J. M. Epstein and R. L. Axtell, Growing Artificial Societies (The MIT Press, Washington, 
D.C., 1996). 

[4] V. Grimm and S. F. Railsback, Individual-based Modeling and Ecology (Princeton, NJ: Uni- 
versity Press, 2005). 

[5] R. L. Goldstone and M. A. Janssen, TRENDS in Cognitive Sciences 9, 424 (2005). 
[6] E. Bonabeau, PNAS 99, 7280 (2002). 

20 



['. 

[8 

[9 
[10 

[11 
[12 

[13 

[14 

[15 

[16 

[it; 

[18 
[19 



R. Axelrod, J. Conflict Resolut. 41, 203 (1997). 

R. Axelrod, The Complexity of Cooperation (Princeton, New Jersey, 1997). 

C. Castellano, M. Marsili, and A. Vespignani, Phys. Rev. Lett. 85, 3536 (2000). 

K. Klemm, V. M. Egufluz, R. Toral, and M. S. Miguel, Phys. Rev. E 67, 026120 (2003). 

K. Klemm, V. M. Egufluz, R. Toral, and M. S. Miguel, Phys. Rev. E 67, 045101(R) (2003). 



L. D. Sanctis and T. Galla. la?Xiv:0707.34 28vl (2007). 

A. Flache and M. W. Macy, |arXiv:physics/0604201| yl (2006). 

D. Jacobmeier, Int. Journ. Modern Phys. C 4, 633 (2005). 

T. Carletti, D. Fanelli, S. Grolli, and A. Guarino, Europh. Lett. 74, 222 (2006). 
F. Vazquez and S. Redner, EPL 78, 18002 (2007). 



N. Boccara, |arXiv:nlin/0611035| yl (2006). 
R. Toral and C. J. Tessone, [aTXiv:physics/0607252| yl (2006). 

D. Centola, J. C. Gonzalez-Avella, V. M. Egufluz, and M. S. Miguel, J. Conflict Resolut. 51, 
2007 (2007). 

[20] C. Gracia-Lazaro, L. F. Lafuerza, L. M. Floria, and Y. Moreno, Physical Review E 80, 046123 
(2009). 

[21] A. Flache and M. W. Macv. larXiv:0808.2710l r2008). 

[22] Y. Shibanai, S. Yasuno, and I. Ishiguro, J. Conflict Resolution 45, 80 (2001). 

[23] J. C. Gonzalez-Avella, M. G. Cosenza, and K. Tucci, Phys. Rev. E 72, 065102(R) (2005). 

[24] J. C. Gonzalez-Avella et al, Phys. Rev. E 73, 046119 (2006). 

[25] J. C. Gonzalez-Avella et al, J. Artif. Soc. Soc. Simul. 10, 9 (2007). 

[26] L. R. Peres and J. F. Fontanari, J. Phys. A: Math. Theor. (2010). 

[27] K. I. Mazzitello, J. Candia, and V. Dossetti, Int. J. Mod. Phys. C 18, 1475 (2007). 

[28] J. Candia and K. I. Mazzitello, J. Stat. Mech. P07007 (2008). 

[29] A. H. Rodriguez, M. del Castillo-Mussot, and G. J. Vazquez, Int. J. Mod. Phys. C 20, 1233 
(2009). 

[30] L. A. Barbosa and J. F. Fontanari, Theor. Biosci. 128, 205 (2009). 



21 



